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I. Introduction 

An edge metric is a metric on the interior of a manifold-with-boundary which is 
singular at the boundary in a manner described by a given fibration of the boundary. 
The related edge differential operators arise in many settings and have been the subject 
of much research; see, for example, |Maj . Our interest in edge metrics arises from the 
observation that they are a robust class of metrics naturally generalizing the product 
of a conformally compact metric with a metric on a compact manifold. The AdS / CFT 
correspondence in physics deals with such product metrics, and the arena of edge 
metrics appears to be a natural setting for the geometric and analytic questions which 
arise. The thesis |Ka] considers a problem concerning eleven-dimensional supergravity 
from this point of view. In this paper we derive a normal form for edge metrics which 
we expect will be useful in further studies. The normal form is the analogue of geodesic 
normal coordinates relative to the boundary at infinity. 

Edge metrics reduce to conformally compact metrics in the special case that the fibers 
of the boundary are points. Asymptotically hyperbolic (AH) metrics are conformally 
compact metrics satisfying a particular scalar normalization at infinity. The normal 
form for AH metrics was derived in |GL] and a different proof was given in |JS2j . This 
normal form has been useful in a number of problems concerning AH metrics. The 
existence statement is that if g is AH on X, then there is a diffeomorphism ip from a 
neighborhood of {0} x dX in [0, oo) x dX to a neighborhood of dX in X such that 
ip\dx = Id and 

/ 1 1 \ / * dx + h x 
1-1 ^9 = 2 ' 

ar 

where x is the coordinate in [0, oo) and h x is a 1-parameter family of metrics on dX. 
The normal form for g is not unique and this is a crucial point. There is a conformal 
class of metrics on dX, called the conformal infinity of g, and the normal forms for g are 
parametrized precisely by the representative metrics in the conformal infinity: h is the 
given conformal representative. Conformal rescalings on the boundary thus correspond 
via the normal form to diffeomorphism changes on the interior. One application of 
the AH normal form is to the renormalized volume of an AH Einstein metric (see 
[G]). The renormalized volume is defined in terms of the exhaustion determined by a 
defining function x in a normal form, and its invariance or noninvariance (measured by 
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a conformal anomaly) is determined via properties of the diffeomorphism determined 
by a conformal change on the boundary. In [JS2J the normal form (II. ip arose in an 
inverse scattering context, where it was used to normalize the action of diffeomorphisms 
on AH metrics. 

For edge metrics there is an analogue of the AH scalar normalization, which defines 
what we call a normalized edge metric. There is another necessary condition for an 
edge metric to have a normal form which is vacuous in the conformally compact case. 
It is that a 1-form on the boundary fibers constructed out of the metric and a defining 
function be globally smoothly exact; this condition is independent of the choice of 
defining function. We call an edge metric satisfying this condition an exact edge metric. 
An exact edge metric has a reduced conformal infinity analogous to the conformal 
infinity in the conformally compact case. In the edge case this is a conformal class of 
metrics on TX\qx/V, where V is the vertical bundle of the boundary fibration (i.e. 
the tangent bundle to the fibers), and any two metrics in the conformal class differ 
by a positive function which is locally constant on the boundary fibers. Our main 
result, Theorem 12.91 asserts that an exact, normalized edge metric g can be put into 
normal form near the boundary by a diffeomorphism which restricts to the identity on 
the boundary, and the different normal forms for g are parametrized precisely by the 
representatives for the reduced conformal infinity. Thus the result takes the same form 
as the result for AH metrics. 

The proof of the normal form for AH metrics in [GLj is based on the observation that 
if ip puts g into normal form as above, then x := x o satisfies the eikonal equation 



By first solving this equation on X, one can therefore obtain directly the x-component 
of . And once one has x, it is clear from the normal form (11.11) that the full map ip 
can be constructed by following the integral curves of the noncharacteristic vector field 
x~ x X$, where X% is the vector field dual to M with respect to g. Now (II. 2p appears to 
be a singular equation. But by writing x = e^x^, where Xq is a fixed defining function 
and bj a new unknown, then expanding the left-hand side of (II. 2p and moving a term to 
the right-hand side, and finally dividing by xq, it becomes a noncharacteristic first order 
nonlinear pde for oj, which is solvable by the method of characteristics. This method 
of reducing a singular eikonal equation to a noncharacteristic initial value problem 
at infinity and then constructing the diffeomorphism by following integral curves of 
an associated noncharacteristic vector field can be used to derive normal forms in a 
number of other settings. It gives a simple alternate derivation of the normal form 
for scattering metrics proved in [JSlj . and is used in [GSJ to derive a normal form for 
0-metrics. (In the case of scattering metrics, the corresponding eikonal equation is 
|§§|^ = 1. One writes x = x + uxl for a new unknown u and divides the equation 
by xl rather than xq to obtain a noncharacteristic initial value problem for u. The 




dx 2 



1. 
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diffeomorphism is constructed by following the integral curves of x~ 2 X$, where X% is 
the vector field dual to |§ with respect to g.) 

In the edge case we were unable to reduce (11.2j) to a noncharacteristic problem. The 
issue is the following. Equation (ll.2p involves the components of the inverse metric g^ 1 . 
In the cases of conformally compact or scattering or 0-metrics, all components of g~ x 
vanish at the boundary. Because of this, one can divide (II. 2p by the correct power of 
Xq as described above to obtain a noncharacteristic problem, and still have an equation 
with smooth coefficients up to the boundary. But for edge metrics, the components of 
g^ 1 along the fibers do not vanish at dX. The coefficient of the derivative transverse to 
the boundary still vanishes there, and one cannot carry out the division by xq to make 
the problem noncharacteristic. So it seems that one is forced to deal with a singular 
equation. 

Joshi studied the normal form for b- metrics in [J]. A b- metric is an edge metric in 
the opposite extreme case from a conformally compact metric: the case in which there 
is only one boundary fiber, the boundary itself. Although geometrically simpler than 
general edge metrics, ^-metrics already exhibit the essential analytic difficulty of the 
general far as derivation of the normal form. By formal calculations, Joshi 

derived the normal form for exact, normalized b- metrics modulo error terms vanishing 
to infinite order at the boundary using the method of [JST], [JS2], but was unable to 
obtain the result in an open set. In this regard he commented, "it is not clear how to 
proceed". The normal form modulo infinite-order vanishing error terms was derived 
for general exact, normalized edge metrics in [Kaj by formal analysis of (jl.2p . (We 
remark that the literature concerning b- metrics, for instance [Mej and [J], defines an 
exact b- metric to be what we refer to as an exact, normalized b- metric. That is, we 
separate the two separately invariant conditions: exactness and normalization.) 

To solve (11. 2p in the edge case, we derive a result concerning existence and uniqueness 
of certain characteristic nonlinear first-order initial value problems for a real scalar 
unknown. We consider problems of the form 

(1.3) xd x u = F(x,y,u,d y u), u(0,y) = u (y). 

Here (x,y) are coordinates on IR n , n > 1, with x > 0, y G M n_1 , and the unknown 
function cu(x,y) is real-valued. F(x,y,co,q) is a smooth real function of (x,y,u,q) G 
R x ]R n_1 x R x ]R n_1 and the initial value Uq is a smooth function of y. 

Theorem 1.1. Suppose that for ally one has 



(1.4) 



F(0,y,Uo(y),d y u) (y)) = 0, 



F q (0,y,u) (y),d y u (y)) = 



and 



(1.5) 



Fw(0,y,wo(y),dyU (y)) < 1. 



Then there exists a unique smooth solution of (11.31) for sufficiently small x > 0. 
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We prove Theorem 11.11 by an adaptation of the method of characteristics. In the 
noncharacteristic case this proceeds by solving the ordinary differential equations cor- 
responding to a Hamiltonian flow-out in the first jet bundle of the solution. Because 
our initial value problem is characteristic, the Hamiltonian vector field vanishes identi- 
cally on the initial submanifold, so there is no flow-out in the usual sense. Nonetheless, 
we are able to construct a flow-out by considering what we call characteristic integral 
curves emanating from a zero of a vector field. We show these exist under an appro- 
priate hypothesis on the eigenvalues of the linearization of the vector field at the zero 
(Theorem 13. 3j) and the union of the characteristic integral curves of the Hamiltonian 
vector field starting from the initial submanifold gives the 1-jet of the solution. It seems 
likely that Theorem 1 1 . 1 1 and the method of using these characteristic integral curves will 
be useful in other problems. Once we have solved ( ll.2p . the second part of the proof 
of the normal form, constructing the diffeomorphism by flowing along integral curves, 
works just as before: the relevant vector field x~ 1 X$ is smooth and noncharacteristic. 

In §2 we define edge metrics and study the geometric structure they induce at infinity. 
We work with edge metrics of arbitrary signature, under an additional nondegeneracy 
hypothesis which we call horizontal nondegeneracy. This is automatic in the case of 
definite signature. We define the notions of normalization, exactness, and reduced 
conformal infinity referred to above, and show that an arbitrary horizontally nonde- 
generate edge metric invariantly induces a metric on each fiber of the boundary. We 
formulate the normal form condition and state the main result, Theorem 12. 9[ which 
asserts the existence and uniqueness of the normal form. We show that Theorem 12.91 
follows from the solvability of the eikonal equation, and reduce the eikonal equation to 
a problem of the form (jl.3p . In §3 we prove Theorem 1 1.1 1 

Throughout, by smooth we mean infinitely differentiable, and all objects are assumed 
smooth unless explicitly stated otherwise. 

2. Edge Metrics 

Let (X, dX) be a manifold-with-boundary. Suppose that dX is the total space of a 
fibration 

F y dX 

TT 

Y 

with fiber F and base Y . One says that (X, dX) is an edge, or boundary-fibered, 
manifold. A motivating special case is that of a product X = M x F, where (M, DM) 
is a manifold-with-boundary. Then dX = DM x F with projection 7r : dX —y DM = Y. 

Let V = ker7r* C TdX denote the vertical vector bundle over dX. A vector field £ 
on X is said to be an edge vector field if its restriction to dX is tangent to the fibers. 
Equivalently, one requires that £\qx G r(V). 

Near a point of dX one can choose local coordinates (x, y a , z A ) so that x is a defining 
function for dX with x > in X°, y a are coordinates on X whose restrictions to dX 
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are lifts of local coordinates on Y (so the y a are constant on the fibers), and z A are 
coordinates on X whose restrictions to each fiber of dX are local coordinates on the 
fiber. The edge vector fields are then spaxi Coo r X \{xd x , xd y a, 8 z a}. 

There is a vector bundle e TX on X, the edge tangent bundle, which can be charac- 
terized by the requirement that its smooth sections are the edge vector fields. xd x and 
the xd y a define sections of e TX which are nonvanishing on dX, and in local coordinates 
near dX one has e T p X = speai{xd x , xd y a, 8 z a} for p £ X. The dual edge cotangent 
bundle e T*X has fibers spanned by the dual basis — , — , dz A . There is a well-defined 
evaluation map Eval : e TX — > TX with the property that Eval p : e T p X — > T P X is an 
isomorphism for p £ X \ dX, but Ev&l p ( e T p X) = V p for p £ dX. For p £ dX, we define 
the horizontal bundle H p = kerEval p C e T p X, so that H p = span{xd x , xd y c}. 

An edge metric g on an edge manifold X is a smooth nondegenerate section of 
S 2 { e T*X). (We allow metrics of arbitrary signature.) In local coordinates near dX, if 
we write 

(<?oo 9o/3 9ob \ / if 
9a0 §<# 9aB 
9A0 9AP gAB j \ dz B 

then this is the requirement that the ^-matrix be smooth and nondegenerate up to dX. 
An edge metric restricts to a usual metric on X° which is singular at dX with the form 
of the singularity determined by the boundary fibration. In the special case that the 
fiber F is a point, there are no z variables and an edge metric is just a conformally 
compact metric on the manifold- with-boundary (X, OX). In the special case that Y 
is a point, there are no y variables, and in this case edge metrics are called b- metrics 
on the manifold-with-boundary (X,dX). For a product edge manifold X = M x F 
as above, a product edge metric is a metric of the form g = gu + gF, where gu is a 
conformally compact metric on (M, dM) and gp is a metric on F . 

We will say that an edge metric is horizontally nondegenerate ii g\}j v is nondegenerate 
for all p £ dX. Clearly any positive definite edge metric is horizontally nondegenerate. 
A horizontally nondegenerate metric induces a metric on the bundle H C £ TX\q X - 
This induced metric may be reinterpreted as a conformal class of metrics on TX\qx/V 
as follows. Multiplication by a defining function x induces a map T p X — > e T p X for 
p £ dX with kernel V p and range H p . This induces an isomorphism H = TX\qx/V 
dependent on the choice of x only up to scale. Via this isomorphism we can transfer the 
metric g\n to a conformal class of metrics on TX\qx/V , called the conformal infinity 
of g. (This choice of terminology is slightly at odds with usual usage in the conformally 
compact case, where the conformal infinity is typically regarded as a metric on TdX 
rather than on TX\qx- One reason for our choice is that TdX/V has rank for the 
special case of 6-metrics.) Of course one can realize the representatives in the conformal 
infinity directly without recourse to H: the metric (x 2 g)\x° extends smoothly to X as 
a section of S 2 T*X and {x 2 g)\gx annihilates V, so it induces a quadratic form on the 
bundle TX\qx/V over dX which is nondegenerate iff g is horizontally nondegenerate. 
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There is a generalization to edge metrics of the normalization condition that a con- 
formally compact metric be asymptotically hyperbolic. The section — of e T*X re- 
stricts to a section of H* which is independent of x. This is because if x = ax with 
< a G C°°(X), then # = — + — , and the restriction of — to H vanishes. If q 
is horizontally nondegenerate, then we can consider the length squared of —\h with 
respect to g\u, and this will be an an invariant of g. 

Definition 2.1. An edge metric g is said to be normalized if g is horizontally nonde- 
generate and if 

dx 2 
— „ =1 on OX. 

Remark 2.2. In the case of indefinite signature, one could equally well consider the 
condition l^fl^L = — 1- Our treatment applies to this case upon replacing g by —g. 

We can also consider the length squared of — with respect to g on all of e TX. In 
general this will depend on x since — \qx is a nontrivial section of e T*X\gx if a varies 
along the fibers. 

Definition 2.3. Let g be an edge metric and x a defining function for dX. x is said 
to be q -normalized if \ — \ 2 = 1 on dX. 

We make these invariant conditions explicit in local coordinates. If g is written as 

( 12. ip . then g\n is represented by ( ^ 00 E ^ ) . Horizontal nondegeneracy of g is the 



9a0 9a, 

requirement that this quadratic form be nondegenerate at dX. In this case, the dual 
metric is given by the inverse matrix, so we write 

#00 9o, \ 1 _ ( C * 

9a0 9a/3 J V * * 

and C = I — I „| 2 , . Thus q is normalized means exactly that C = 1, and this condition 

I x \H\g\ B y J 

is independent of the choice of all the coordinates. 

On the other hand, — is a dual basis vector in the full frame. We write 

' X 

900 9ob V 1 / 9°° 9°" 9° B 

9a, 9aB \ = \ r° 9^ 9^ 
n 7, \ \ nAO z,AB 

9a, gAB j \ g 9 9 

and then ^ 00 = | ^ \ 2 . So x is ^-normalized means g 00 = 1. This condition is independent 
of the choice of y a , z A , but in general does depend on the choice of x. 

A horizontally nondegenerate edge metric invariantly induces a pseudo-Riemannian 
metric on the fibers of dX. Let g be an edge metric. The induced dual metric on e T*X 
is a section g _1 of S 2 ( e TX). Thus Eval(g -1 ) is a smooth section of S 2 TX. On dX, this 
section of S 2 TX degenerates: its restriction to dX is a smooth section of S 2 V C S 2 TX. 
Elementary linear algebra (most easily carried out in terms of the explicit formulation 
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of these conditions below) shows that the condition that g is horizontally nondegenerate 
is equivalent to the condition that Eva^g" 1 ) is a nondegenerate section of S 2 V. So if g 
is horizontally nondegenerate, Eva^g™ 1 ) defines a metric on V*. Its dual is a metric on 
V, or equivalently a pseudo-Riemannian metric on each fiber of dX. We denote this 
induced metric on the fibers by gp. 
Concretely: 

g 00 g^ g 0B \ ( xd x 
(2.2) g- 1 = ( xd x xd ya 8 z a ) I g a0 g afi g aB xd yl 



gAO - g A P gAB \ d z 



so Eval(g 1 )\qx = g AB d z Ad z B . Nondegeneracy of ( ^ 00 ^ 0/3 ) is equivalent to nonde- 

\ 9a0 9a/3 J 

generacy of g AB '. The induced metric on the fibers is [gp)ABdz A dz B , where (gp)AB = 
(g This metric gp is independent of the choice of all coordinates. 

Next we introduce the notion of an exact edge metric. Let g be an edge metric 
and x a defining function. Now — is a smooth section of e T*X. Let X x be the edge 
vector field dual to — with respect to g. EvaX(X x )\ax is then a section of V. If g 
is horizontally nondegenerate, we can define a section a x of V* to be its dual with 
respect to gp. Thus we have associated to each defining function x a 1-form a x on the 



a. 



fibers of dX. If x = ax is another defining function, then ^ = ^£ + ^ = ^ + rflog 

° ' x x a x 

Following through the definition shows that the corresponding 1-forms are related by 
otax = «i + dy log a, where dyloga = d\oga\v- Thus a x changes by an exact form 
under change of defining function. 

Definition 2.4. An edge metric g is said to be exact if g is horizontally nondegenerate 
and if for each defining function x, there is / G C°°(dX) so that a x = dyf. 

The above reasoning shows that if this holds for one x, it holds for all x. If g is exact, 
then by correct choice of a we can find x so that a x = 0. 

Definition 2.5. If g is an exact edge metric, a defining function x is said to be g-related 
if a x = 0. 



If x is g-related, then another defining function x = ax is also g-related if and only if 
dya = 0, i.e. a\gx is locally constant on the fibers. So g-related defining functions are 
determined precisely up to multiplication by a positive function whose restriction to 
dX is locally constant on the fibers. 

Clearly x is g-related if and only if Eval(X x )\sx = 0. Since g" 1 is given by (12.21) . we 
deduce that x is g-related if and only if g 0B \ox = 0, and this condition is independent 
of the choice of y a and z A . 

Recall that a horizontally nondegenerate edge metric g induces a conformal class of 
metrics on TX\q X /V with representatives (x 2 g)\sx- If 9 is exact, we can restrict to 
representatives of the conformal class which arise from g-related defining functions x. 
We will call this the reduced conformal infinity of g. Metrics in the reduced conformal 
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infinity are determined up to rescaling by functions which are locally constant on the 
fibers of dX. Choosing a representative metric in the reduced conformal infinity is 
entirely equivalent to choosing a g-related defining function x to first order at dX. 
The following lemma will be useful in the sequel. 

Lemma 2.6. Let g be an exact edge metric and x a g -related defining function. Then 
\^r\ H \ g \ H — \^~\g on dX. In particular, if g is exact and normalized, then every g- 
related defining function is g -normalized. 

Proof. Choose local coordinates (x, y a , z A ), taking x to be the given ^-related defining 
function. Recall that if 

' T U 

V w 

is an invertible matrix in block form with T and W square and T invertible, then its 
inverse can be written 

T- 1 + T^US^VT- 1 -T^US- 1 
-S- l VT- 1 S- 1 



M 



M- 1 = 



where S = W — VT 1 U necessarily is invertible. Apply this with 

(#00 90f3 9ob 
9a0 9af3 9aB 
9A0 (jA/3 ()AB 

\ 9a0 9af3 J \ 9aB J 

all evaluated at dX. The hypothesis that x is g-related says exactly that the first row 
of T^US^ 1 vanishes. Hence the first row of T^US^VT^ 1 also vanishes. Hence the 
first row of T _1 + T^US^VT^ 1 agrees with the first row of T _1 . In particular, their 
00 components agree, which is the desired conclusion. □ 



Remark 2.7. A simpler proof can be given if one assumes that \ — \ 2 ^ on dX. This 



X iff 



is of course automatic for g positive definite. Under this hypothesis one can make a 
change of the y- variables y a — >■ y a + \ a x, with \ a \ox chosen to make g 0a \ox = 0. In 
the new coordinates one has g 0a = and g 0A = on dX, and the conclusion is clear. 

Next we formulate the normal form condition. If X is an edge manifold, then dX is 
the total space of a fibration. Consider [0, oo) x dX as a manifold-with-boundary, with 
boundary {0} x dX = dX. The given fibration of dX induces a natural edge manifold 
structure on [0, oo) x dX. The coordinate x of the first factor is a canonical defining 
function on [0, oo) x dX. 

Definition 2.8. An edge metric g on a neighborhood U of {0} x dX in [0, oo) x dX is 

in normal form if g — + k, where A; is a smooth section of e T*U satisfying xd x jk = 
everywhere. 
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This is equivalent to requiring that g have the form 



9=(~ dJ T dzA )\ 9a0 9*3 



I 1 



\ g Af 3 9ab 



) 




) 



and k is given by 



k = ( 4£ rf^ ) 



#0$ S'aB 
#A/3 9AB 



) 




) 



Observe that g is exact and normalized. Also x is ^-normalized and g-related. 

The main theorem asserts that any exact, normalized edge metric g can be put into 
normal form, and the normal forms for g are parametrized by the g-related defining 
functions to first order, or equivalently by the representatives for the reduced conformal 
infinity. 

Theorem 2.9. Let X be an edge manifold and g an exact, normalized edge metric. 
If Xq is a g-related defining function, then there is a unique diffeomorphism ip from a 
neighborhood of {0} x OX in [0, oo) x dX to a neighborhood of dX in X , such that 
^\dx = Id, ip*g is in normal form, and ip*xo = x + 0(x 2 ). 

The main step in the proof is to solve the eikonal equation: 

Proposition 2.10. Let X be an edge manifold and g an exact, normalized edge metric. 
If xq is a g-related defining function, then in a neighborhood of dX there is a g-related 
defining function x, uniquely determined by the conditions 



Theorem 12.91 follows from Proposition 12.101 by the usual argument of flowing along 
integral curves: 

Proof of Theorem \2.9\ Let x be as in Proposition 12.101 Recall that X% is the edge 
vector field dual to dx/x with respect to g, and Eval(X^) = on dX since x is g- 
related. Consequently iV := x~ l Eval(X^) is a smooth vector field up to dX, and 
Nx = = 1. In particular, iV is transverse to dX. For x > and p G dX, define 
i/j(x,p) to be the point obtained by following the integral curve of iV emanating from p 
for x units of time. Since Nx = 1, we have ip*x = x, and iV is orthogonal to the level 
sets of x since X$ is dual to dx/x. Thus ip*g has the desired form. □ 

We conclude this section by reducing Proposition 12.101 to the solution of a singular 
initial value problem of the form considered in Theorem 11.11 It suffices to prove Propo- 
sition |2HH] locally in a neighborhood of a boundary point, since the uniqueness implies 
that the local solutions will piece together to form a global solution. Relabel xq as x 



dx 2 




1 



X g 
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i:=i- 


dx 




X 



and write x = e^x. Our new unknown is uj, with boundary condition uj = at x = 0. 
Now ^ = — + dou, so the equation |45| 2 = 1 becomes 

x x ' ^ I £ I g 



(2.3) 2X^+^0; 



where we now neglect the distinction between X x and EvaA(X x ). Lemma [2TB1 shows that 
x is g-normalized, so the right-hand side vanishes at dX. Work in local coordinates 
(x,y a ,z A ) as above. The left-hand side is a quadratic polynomial in xd x uj, xd y acu, and 
8 z auj with no constant term and with coefficients smooth up to the boundary. It follows 
that (12.31) can be written as 

(2.4) Q(x, y, z, xd x uj, d y uj, d z uj) = f(x, y, z), 

where Q is a quadratic polynomial in (xd x uj, d y uj, O z uj) with no constant term and with 
coefficients depending on (x, y, z) which are smooth up to x = 0, and / is smooth with 
f(0,y,z) = 0. (We have absorbed the x multiplying d y uj into the coefficients.) We 
make the following observations about Q. First, the coefficient of the linear term xd x uj 
is nonzero at x — 0, since X r x = x\ — \ . Second, the coefficients of the linear terms 

I X I g 

dyU and d z uj vanish at x = 0, since x is g-related so that X x = at x = 0. Third, all 
of the arguments (xd x u , d y u , d z u) themselves vanish at x = when evaluated on any 
function u satisfying the initial condition u = at x = 0. In particular, the partial 
derivative of the quadratic terms of Q with respect to any of xd x co, d y u, d z u vanishes 
at x = when evaluated on the initial data. 

The implicit function theorem (or the quadratic formula) implies that in a neighbor- 
hood of (x, y, z, xd x u, d y u, d z u) = (0, y, z, 0, 0, 0), (12.4R may be solved for xd x u. So it 
may be written in the form 

(2.5) xd x u = F(x,y,z,d y u,d z u), 

where F is a smooth function of its arguments satisfying F(0, y, z, 0, 0) = 0. More- 
over, the observations above show that Fq ^(0, y, z, 0, 0) = and Fq zUJ (0, y, z, 0, 0) = 0. 
Equation (12.51) with initial condition u = at x = is of the form (11.31) . where y in ( 11. 3D 
plays the role of (y,z) in ( I2.5p . Condition (ll.5p holds since F in (I2.5P is independent 
of uj, so that F^ = 0. Thus Proposition 12. 101 follows from Theorem 11.11 



3. Singular Initial Value Problems 

In this section we prove Theorem 11.11 First observe that the conclusion in Theo- 
rem [1J] fails without the hypothesis F u (0, y, Wo(y), d y cJo(y)) < 1. For instance, the 
equation xd x u = u has infinitely many smooth solutions co = cx satisfying cu(0) = 0, 
and the equation xd x u = u + x has no smooth solutions (the general solution is 
co = cx + xlogx). Also note that if < F u (0,y,cj o (y),dyCJ (y)) < 1, then the smooth 
solution need not be the only continuous solution. For example, if < a < 1, then 
uj = cx a solves xd x uj = au with Co>(0) = for any cel. In this case the unique smooth 
solution is uj = 0. 
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We first use a standard reduction technique via Taylor expansion to reduce the 
equation to a simpler form. In the following we denote F^ (y) = F u (0, y, u)o(y), d y uo(y)) 
and similarly for other derivatives of F evaluated on the initial data. 

Observe first that differentiating (11 .3p with respect to x at x = and solving for u x 
shows that if a; is a smooth solution, then 

F^\y) 

u x (0,y) = x :=ui(y). 
1 - F^ >{y) 

We can write 

u(x, y) = uj (y) + x{ui(y) + u(x, y)) 
for a smooth function u(x,y), and regard u as the new unknown. 

Proposition 3.1. In terms of u, ( 11.31) becomes 

(3.1) xd x u = (F^ 0) (y) - l)u + xG(x,y,u,d y u), u(0,y) = Q, 
where G is a smooth function of (x, y, u, q). 

Proof. It is clear from the discussion above that the initial condition on u is w(0, y) = 
0. Set qo{y) = d y u (y). The second order Taylor expansion of F(x,y,ou,q) about 
(0,y,u (y),qo(y)) takes the form 

F(x, y, u, q) = F$>\y)x + F^(y)(oj - cu (y)) + Q(x, u - u (y), q - q (y)), 

where Q is a homogeneous quadratic polynomial of its arguments with coefficients 
which are smooth functions of (x, y, u, q). We have 

d x u — coi + (xd x + l)u, dyU — qo(y) = x(d y ui + d y u). 

Substituting and then dividing by x shows that (I1.3P becomes 

ui + (xd x + l)u = Fi 0) + i^ 0) (wi + u) + xG(x, y, u, d y u) 

for a smooth function G. The definition of u% shows that u>\ = F x ^ + Ui, so this 
reduces to ( 13. lft . □ 

Proposition 13. II implies that Theorem 11.11 follows from the following special case. 

Proposition 3.2. Let b(y) and G(x, y, u, q) be smooth and suppose b(y) < 0. Then the 
IVP 

(3.2) xd x u = b(y)u + xG(x,y,u,d y u), u(0,y) = 
has a unique smooth solution for sufficiently small x > 0. 

We prove Proposition 13.21 by an adaptation of the method of characteristics. The 
main tool is a result asserting the existence and uniqueness of smooth "characteristic 
integral curves" of time-dependent vector fields vanishing at an initial point. 

Let M be a smooth manifold and po G M. Suppose that V(t,p) is a smooth time- 
dependent vector field defined for t near and p in a neighborhood of po, such that 
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V(0,po) = 0. By a characteristic integral curve for V at p we mean a curve 7 : [0, e) — > 
M for some e > such that 

(3.3) t j t ~f(t) = V(t, 1 (t)), 7 (0) =p . 

Recall that the linearization of a vector field at a zero is the endomorphism DV of 
T po M such that V(p) = DV(p — p ) to first order at p . For a time- dependent vector 
field this refers to the linearization of the vector field in the space variables with t fixed. 

Theorem 3.3. Let V(t,p) be a smooth time- dependent vector field such that V(0,po) = 
0. Suppose that all eigenvalues A of DV(0,po) satisfy Re A < 1. Then on a sufficiently 
small time interval there exists a unique smooth characteristic integral curve for V at 
Pq. This characteristic integral curve depends smoothly on variations of the initial point 
Po for which the conditions V(0,po) = and Re A < 1 continue to hold. 

Observe that the case M = R of Theorem 13.31 coincides precisely with the special 
case n = 1 of Theorem II. 1[ upon relabeling t as x, 7 as u>, and V as F. In particular, 
the examples above show that existence and uniqueness of smooth solutions can fail if 
A = 1, and there may be continuous solutions which are not smooth if < A < 1. 

The first step in the proof of Theorem l3.3l is to perform a Taylor expansion analogous 
to the one made above for the pde. Work in local coordinates on M and let Xq denote 
the coordinates of p . Differentiating (13. 3p with respect to t at t = and solving for 
7'(0) gives 

Y(0) = [/ - DV(0, z o )rV t (0, aro) := 7i- 

Write 

(3.4) 7 (t) = x + * (71 + • 

Upon Taylor expanding V(t, x) about (0,x ), substituting (13.4)) . and simplfying as in 
the proof of Proposition 13.11 one finds that when written in terms of a, (13.31) takes the 
form 

(3.5) to-' + Aa = tG(t,a), <r(0) = 0. 

Here A = I — DV(0, Xq) has the property that all of its eigenvalues have positive real 
part. G is smooth, and A and G depend smoothly on the initial point xq. Initial value 
problems of the form (13. 5p are studied in Chapter 5 of |Kij . The results formulated 
there assume that A is independent of the parameters, but the same arguments apply 
to our situation. We briefly outline a proof that ( I3.5P has a unique smooth solution 
varying smoothly with the parameters xq if the eigenvalues of A have positive real part. 
Theorem 13.31 is then a consequence by the reduction above. 

The problem (13.51) can be reformulated as the integral equation 

(3.6) ait) = (Ta)(t) := t [ s A G( y st,a( y st))ds. 

Jo 

The hypothesis that the eigenvalues of A have positive real part implies that the oper- 
ators s A are uniformly bounded for < s < 1. A standard contraction mapping/fixed 
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point argument proves the existence and uniqueness of a continuous solution. To es- 
tablish smoothness in t, rewrite (I3.6j) as 

a(t) =T A [ s A G{s,a{s))ds. 
Jo 

This shows that a is C 1 for t > 0. Differentiate in t and change variables back to obtain 

(3.7) a\t) = G(t, c(t)) -A I s A G(st, a(st))ds. 

Jo 

Thus a is C 1 up to t = 0. Now successively differentiating (13. 7p shows that a is 
C°°. Smoothness of a with respect to the parameters is a consequence of the implicit 
function theorem applied to the equation a — Ta = 0. 

Proof of Proposition \3.2[ We construct a singular version of a Hamiltonian flow-out in 
the first jet bundle of the solution u. The argument follows the usual reasoning for the 
non-characteristic case, substituting Theorem 13.31 in an appropriate parameterization 
for the existence and uniqueness of integral curves of the Hamiltonian vector field. 

Let J denote the first jet bundle of a scalar function u on R n , with coordinates 
(x,y,u,p,q) where p is the variable dual to x, and projection it : J — )■ W 1 given by 
ir(x, y, u,p, q) = (x, y). Set x = (x, y) and p = (p, q). The 1-jet of a function u on W 1 is 
the section of J given by x — > (x, w(x), <iw(x)). We denote its image {(x, w(x), du(x))} 
by S u ; this is a submanifold of J of dimension n. The tautological contact form is 
9 = du — pjdx\ S u is a Legendrian submanfold relative to 8; i.e. the pullback of 8 to 
S u vanishes. 

Recall that if -ff(x, u, p) is a smooth real function on J ', the associated Hamiltonian 
vector field is 

£,h = H Pi d^i + ~piH Pi d u — (H x i + piH u )d Pi . 
It is uniquely determined by the conditions 

(3.8) J d0 = dH mod 8 0(£ H ) = 
and satisfies 

(3.9) £ H Ad8 = dH- H u 8, £ H H = 0. 

If u is a solution of H (x, u, du) = 0, then is tangent to S u at all points of S u . 
Take to be the Hamiltonian corresponding to (13.21) . i.e. 

H(x, y, u, p, q)=xp- b(y)u - xG(x, y, u, q). 

Differentiating (13. 2p at x = shows that a solution u must satisfy 

f*m\ /» ^ G(0,?/,0,0) 

(3.10) 9«tt(0,y)= : =PoW- 

Define a smooth submanifold X C J of dimension n — 1 by 

X={(0, 2 /,0,p (y),0)}. 
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Hamilton's equations for the integral curves of take the form 

dx 



x 



— = -xG n 



ds 
dy 

ds 
du 

ds 
dp 

ds 
dq 

ds 



x(p - qiG qi 

p( b (y) - i) 



G 



x(G x +pG u ) 
^ = b y (y)u + b(y)q + x(G y + qG u ). 



Observe that vanishes identically on X. So all integral curves of beginning 
on X are constant; there is no Hamiltonian flow-out in the usual sense. Instead we 
consider characteristic integral curves of beginning on X. is time-independent 
and constants are also characteristic integral curves. But the characteristic integral 
curves are not unique: D^h on X has A = 1 as an eigenvalue arising from the first 
equation in the system above. By using x as the parameter, we will obtain unique 
nonconstant characteristic integral curves of emanating from X whose union will 
form the submanifold S u giving the solution u. 

Use i as a parameter for the characteristic integral curves. The first equation above 
gives d/ds = xd/dx. Substituting in the remaining equations gives the system 

dy 



x— = —xG„ 
dx 

du 



x 



(3.11) 



x 



x(p-qiG qi ) 

p(b(y)-l) + G + x(G x + P G u ) 



dx 
dp 

dx 
dq 

x— = b y (y)u + b(y)q + x(G y + qG u ). 

This has the form (13. 3p . where x plays the role of t. Choose yo G 
initial conditions 



su- i 



and impose 



(3.12) 



1/(0) =3/0, «(0)=0, p(0)=po(2/o), 9(0) =0. 



The linearization of the right-hand side of (13. lip evaluated at x 
initial conditions for (y,u,p,q) is 



and at the given 



D 



/0 




\ 



Kvo) - 1 * 

b(y )lj 
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where the blocks have sizes n — 1, 1, l,n — 1. Here * denotes a quantity whose value 
will be irrelevant and I denotes the {n — 1) X (n — 1) identity matrix. The eigen- 
values of D are with multiplicity n, b(yo) — 1 with multiplicity 1, and b(y ) with 
multiplicity n — 1. These are all real and less than 1, so Theorem 13.31 implies that 
there is a unique smooth solution (y(x,yo),u(x,yo),p(x,yo),q(x,yo)) of (13.111) . (13.121) 
for sufficiently small x varying smoothly with y^. 
Define a map $ into J by 

$(x,yo) = (x,y(x,y ),u(x,y ),p(x,y ),q(x,y )). 

Since 2/(0, 2/0) — Vo, it follows that $ is a diffeomorphism from a neighborhood of 
{x = 0} to a submanifold T C J of dimension n. We claim that H = on J 7 . Since 
£h% = x and yo is constant on the solution curves, it follows that — X 9 X - Since 

£hH = 0, we have xd x (§*H) = 0. Since H = on X, one concludes that iJ = on J 7 
as claimed. 

We now prove existence in Proposition 13.21 The projection 71 : J — > M. n restricts to 
a diffeomorphism from J 7 to a neighborhood of {x = 0} (possibly after shrinking J 7 ). 
Therefore on J 7 we can regard u,p, q as functions of (x, y). In particular this defines a 
smooth function u(x,y). We claim that 

(3.13) p(x, y) = d x u(x, y), q^x, y) = d y m(x, y). 

This is equivalent to saying that J 7 = S u . Existence in Proposition 13.21 follows imme- 
diately, as then the equation H = on T together with the initial condition become 
the statement that u satisfies (13.21) . 

Since 9 = du — pdx — qidy 1 , in order to prove (13.131) it suffices to show that the 
pullback of 9 to J 7 vanishes. Recalling (13.81) . (13.91) . observe that 

C iH 9 = i H jd9 + d(0(fn )) = dH - H u 9 } 

where £ denotes the Lie derivative. For the pullback to J 7 we therefore obtain C^ H 9 = 
-H u 9. Pulling back by $ gives C xdx $*9 = - ($*H U ) If we write $*9 = 9 dx + 
9id(y ) 1 , then this becomes 

(3.14) xd x 9 = (-$*H U - 1) 9 , xd x 9 t = t . 

For each y , these are scalar ode's of the form xd x v = 0(x)v, where v = 9 or 9i and 
f3 = —§*H U — 1 or —§*H U . Since H u = —b(y) at x = and b < 0, we have (3 < near 
x = 0. The general solution is 

r (3{s)ds m 

v = cexp / = cx^ y > li{x) } 

J s 

where c G K and h is a nonvanishing smooth function. Since /3(0) < and 9 is smooth, 
we must have c = 0, so we obtain 9q = 9 { = as desired. (Alternately, the vanishing 
of 9 and 9i follows from uniqueness in Theorem 13.31 applied to (I3.14p . ) 

Finally we prove uniqueness in Proposition 13.21 We show that if u is any smooth 
solution of (13.21) . then S u = T near X. We have already observed that (^(O, y) is given 
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by frXTO]) . so that S u n {x = 0} = X = J 7 n {x = 0}. The system 

for unknown y(x) has a unique smooth solution (either by Theorem 13. 3l or by cancelling 
x and quoting usual ode theory). As y varies, the corresponding curves (x,y(x)) fill 
out a neighborhood of {x = 0} in IR n fl {x > 0}. Therefore near X, S u is the union over 
y of the lifts 

(3.15) x h> (x, y(x),u(x, y{x)), d x u(x, y{x)), d y u(x, y{x))). 

The curves 

x i y (y(x),u(x,y(x)),d x u(x,y(x)),d y u(x,y(x))) 
solve (13.111) since £ H is everywhere tangent to S u . Since T was defined to be the union 
of all curves (I3.15P corresponding to solutions of (13.111) , it follows that S u = T near X. 
□ 
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